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■ Abstract 

O 
O 

I The trace over the degrees of freedom located in a subset of the space trans- 
om' forms the vacuum state into a mixed density matrix with non zero entropy. 

■ This geometric entropy is beheved to be deeply related to the entropy of black 
. holes. Indeed, previous calculations in the context of quantum field theory, 

where the result is actually ultraviolet divergent, have shown that the geomet- 

^ I ric entropy is proportional to the area for a very special type of subsets. In 

J> ' this work we show that the area law follows in general from simple consider- 
ations based on quantum mechanics and relativity. An essential ingredient in 

CN ■ our approach is the strong subadditive property of the quantum mechanical 

^ i entropy. 
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I. INTRODUCTION 



Oh. a chain of arguments using thermodynamics, classical general relativity and quantum 

r-| ! field theory in curved spaces shows that the black hole has an associated entropy proportional 
to the horizon area. This is arguably one of the most important clues to the understanding 
of the role of gravity at the quantum level. However, these calculations give no precise 
identification of the degrees of freedom responsible for the entropy, nor an explanation of 
its relation with area. There are several approaches in this regard, involving either specific 
models of quantum gravity or effective ideas [1]. 

A very appealing candidate in this last sense is given by the entropy of entanglement 
[2,3]. On the black hole spacetime we can write the Hilbert space of a quantum field as a 
tensor product of two Hilbert spaces, which correspond to the asymptotic data at infinity 
and at the horizon respectively. The quantum state relevant for asymptotic observers is 
obtained by tracing the complete quantum state density matrix over the invisible 

degrees of freedom on the horizon. But contains correlations between the field modes 
inside and outside the black hole, or, in other words, it is entangled with respect to the 
above tensor product. Thus, the partial trace leads to a mixed density matrix with non zero 
entropy. 

For large black holes the near horizon geometry is approximately flat. This suggests 
that the relevant entropy should also occur in flat space [4]. To mimic the black hole in 
Minkowski space take a spatial hyperplane and a subset X on it. This divides the space in 
two, the interior and exterior regions with respect to X. The trace of the vacuum state over 
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the interior (or exterior) degrees of freedom is a mixed density matrix. Its entropy, being 
a function of the subset X, was called geometric entropy [4]. This new flat space quantity, 
suggested by black hole physics, is interesting in its own right, and provides a simpler context 
to contrast ideas. 

Several authors have done explicit calculations of the geometric entropy in quantum field 
theory (QFT). The methods ranged from numerical work [3,5], to analytical calculations 
based on the replica trick, the heat kernel, and conformal symmetry in two dimensional 
spacetime [4,6-10]. In all cases, the reassuring result is an entropy proportional to the area. 
However, most of the work has focused on free fields and a very restricted class of subsets. 
The application of the analytical methods was constrained to the infinite half space since 
they use in an essential way the fact that density matrix corresponding to the Rindler wedge 
is a thermal one [10]. The numerical methods encounter practical problems for non spherical 
subsets (sec however the recent work on annular and flower-like geometries [11,12]). Also, 
the physical principles behind the proportionality between area and entropy remain obscure 
beyond the explicit calculations, that involve dealing with non renormalizable divergences. 

In this work we show that for a wide class of subsets the area law for the geometric entropy 
follows from very general and model independent ideas based on quantum mechanics and 
relativity. 

We assume, as it has been done explicitly or implicitly in all the previous work on the 
subject, that it is possible to divide a Cauchy surface into pieces and trace the vacuum state 
over the degrees of freedom in one of them to obtain a density matrix. This is the single 
basic condition allowing to define a geometric entropy. However, this is a delicate point, 
related to the divergences appearing in the QFT calculations. We will only say a bit more 
about it in Section HI. 

The causal evolution is unitary. Accordingly, the geometric entropy corresponding to 
two space-like subsets coincide if they give place to the same causal domain of dependence, 
or, what is the same, they can be extended to a global Cauchy surface by union with the 
same space-like subset. We also make use of the Poincarc symmetry of the vacuum state. 

It is essential to the arguments below the mathematical fact that the quantum mechanical 
entropy is positive and strongly subadditive. The strong subadditive property was intro- 
duced in [13,14], in the context of the investigations around an old problem in statistical 
mechanics, namely, the proof of the existence of the mean entropy (the limit of entropy over 
volume for large sets). Here we apply these ideas to the relativistic situation. Note that for 
a general quantum statistical mechanical system the entropy in a given volume has to be 
defined in the same way as the geometric entropy is defined above, that is, by tracing over 
degrees of freedom in the appropriate domain. In this sense, the geometric entropy is just 
the name for the more familiar concept of the entropy contained in a volume in the special 
case where the state is the relativistic vacuum. 

The cited investigations on the Euclidean symmetric situation where mainly focused 
on the extensive properties of the entropy on large domains. Here the mean entropy is 
zero. Basically, this is so because we can construct Cauchy surfaces with arbitrarily small 
volume by approaching the null surfaces, and use the boost symmetries and the subadittive 
property to put an upper bound to the entropy which grows like the area. Moreover, the 
strong subadditivity allows to completely determine the functional form of the geometric 
entropy. Curiously, a very relevant piece of the information it provides, that gives the clue 
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to show that the geometric entropy is not only bounded by the area but hnearly growing 
with it, shows up ahcady when combined only with translation invariance. However, the 
crucial limit for the geometric entropy is not the one of infinite volume but that of the very 
flat type of sets. 

Summarizing, we show that the combination of strong subadditivity, Poincare symmetry 
and caTisality constrains the geometric entropy to be proportional to the bounding area, plus 
a constant term. 

The work is organized as follows. In Section II we review some properties of the quan- 
tum entropy. In Section III we state the problem of the geometric entropy in Euclidean 
and relativistic spaces. In Section IV we review some work on the entropy for Euclidean 
symmetric states and develop some new results on this subject, related to the entropy of 
very fiat sets. The main result of the paper is the theorem in Section V which states that 
geometric entropy is proportional to the area. Finally we present the conclusions and discuss 
some open perspectives. 



In this Section we review some inequalities satisfied by the quantum entropy. The proofs 
of the statements, the original references and further details can be found in the reviews 



In quantum mechanics a physical state is described by a density matrix p, which is a 
self-adjoint and positive operator with unit trace, Trp — 1, defined on a Hilbert space Ti. 
The entropy S is function of the quantum state given by the expression 



The quantum entropy is always positive, 5" > 0, being zero if and only if the state is pure, 
that is, the density matrix is a one dimensional projector p = \^){^\. The entropy can be 
infinite for infinite dimensional Hilbert spaces, while on a space with finite dimension d it is 
bounded above by log d. The function S satisfies 



Ai5(pi) + X^S{p2) - Ai log(Ai) - A2 log(A2) > ^(Aipi + A2P2) > Ai5(pi) + As^lps) , (2) 



where Ai and A2 are any positive numbers such that Ai -|- A2 = 1. The second inequality 
means that the effect of mixing always increases the entropy. Moreover, any state can be 
decomposed as a mixing of pure states. 

In the following we will be mainly interested in the properties of the entropy related to 
spaces constructed out of the tensor product of two or more Hilbert spaces. Let the Hilbert 
space 7i = Til 07^2 and p = Pi^i2,i\i'^ be the density matrix on H.. By tracing over Tii we 
can form the density matrix p2 = (^2)121^ — Phi2,w'2 The density matrix pi is 

defined similarly by tracing over 7-^2 • 

For a pure density matrix p the following duality property for the entropy of the subsys- 
tems holds 



II. PROPERTIES OF THE QUANTUM ENTROPY 



[15,16]. 



S = —Trp log p . 



(1) 



S{p,) = S{p2) . 



(3) 
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This follows from the fact that for p — |^)(^| we have the equality Tr{p'^) = Tr^p^) for any 
power n, and consequently the non zero eigenvalues and their multiplicities coincide for pi 
and p2- 

It is easy to give examples where p is pure, so that S{p) = 0, while pi is in a mixed state 
with S{pi) > 0. Thus, the quantum entropy can not be generically increasing (monotonicity) 
with the size (the number of degrees of freedom) of the subsystem (nor with the size of the 
subsystem that is traced over). This is somehow against the intuition coming from the 
extensive entropy of gases. 

A natural question is what type of density matrices pi can arise from tracing a pure 
density matrix p over a subsystem. The answer is that any density matrix can be obtained 
in this way. A density matrix pi on a space H can always be realized as the partial trace of 
a pure density matrix p onTi^Ti. 

The entropy for a system composed by two independent subsystems is the sum of the 
subsystem entropies. This is at the physical basis of the explicit form of S given by equation 
(1). The density matrix for the composed system is given by the tensor product p — pi® P2 
in this case, and we have the equation 

S{p) = S{p{) + 5(p2) . (4) 

In the general case where p does not factorize and the subsystems can not be considered 
as statistically independent, this equation can be extended to an inequality expressing the 
subadditive property of the entropy 

S{p) < S{p,) + S{p2) . (5) 

In fact the entropy satisfies a stronger inequality. To introduce it, it is convenient to 
consider a more general case in which the Hilbert space of the system is a tensor product 
of an arbitrary number of factors 7i = 'S'iei'^ij where / is the set of indices labeling the 
different subsystem Hilbert spaces Hi- Let A be any subset of the set of indices /. Define 
the corresponding reduced density matrix on 0jg^ Hi, by tracing over all Hi with i ^ A, 

Call S{A) = S{pa) the corresponding entropy. With this notation the subadditive property 
(5) writes 

S{A) + S{B) > S{AU B) (7) 

for any pair of subsets A and B of /, while the duality relation for a pure total density 
matrix p is simply 

SiA) = Si-A) . (8) 

Here —A is the set complementary to A in /. As mentioned, / can be artificially enlarged 
to a bigger set and the total density matrix can be taken pure in the enlarged space, 
without modifying the value of S{A) for A <0 I. In this case (8) holds in the bigger space /'. 

Applying subadditivity to the complements —A and —B we have 'S'(— ^4) + S{—B) > 
S{—{A (IB)), since the complement operation interchanges union and intersection. Thus, 
purifying the total density matrix the inequality (5) implies a different relation, namely 
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S{A) + S{B)>S{Ar)B). (9) 

The strong subadditivity (SSA) generalizes the relations (7) and (9) in a form that is self-dual 
under taking complements in the pure case. It writes 

S{A) + S{B)>S{AUB) + S{AnB). (10) 

Using the trick of working with an artificially purified density matrix the SSA applied 
to A and —B leads to a different self- dual relation. This is 

S{A) + S{B)>S{A-B) + S{B -A), (11) 

where A — B means Ar\{—B). The inequalities (10) and (11) are equivalent in the pure case 
but they are both valid in general. In lack of a better name, and since it is closely related 
to SSA, we often call SSB the inequality (11). 

Finally, we also mention that for two non intersecting sets A and B the SSA and SSB 
relations lead to the following triangle inequalities for the entropies 

\S{A) - S{B) \<S{AUB)< S{A) + S{B) . (12) 

The second one is of course subadditivity. The inequahties (12) are useful to prove some 
continuity properties of S. For example, if the entropy of B is very small one has that S{A) 
and S{A U B) are very near to each other. 

We mentioned that S is not monotonic in general. By monotonicity here we mean the 
property S{B) < S{A) for any BOA. However, the inequality (11), or the first one from 
(12), can be thought as partial compensations for the lack of monotonicity. 

Inequalities provide less concrete information and are more difficult to handle than equa- 
tions. This is perhaps the reason why the relations (10) and (11) are relatively poorly known 
and have found only few applications. However, in the presence of symmetries these can be 
very useful as will become apparent in the following. 



III. STRONG SUBADDITIVITY AND CONTINUUM QUANTUM SYSTEMS 

Let us consider a continuous quantum system in R*^ described in terms of a Fock space. 
The single particle states are given by elements of TCiilV^) = /^^(R'^), the Hilbert space 
of square integrable functions in R'^, while the n particle space is the product of n copies 
of 7^1 (R"^), 7in(R'*) = 7ii(R'')(8)...(8)7t;i(R'^), where the tensor product is understood as 
symmetrized or antisymmetrized according to the particles being bosons or fermions. The 
Fock space is then given by the direct sum H — 0^ Hn- Similarly, we can define the single 
particle states in a given volume V in R*^ as HiiV) = C'^iV), and its corresponding Fock 
space Ti-lV) = 0^Hn(^)- If the intersection of two sets Vi and V2 has zero measure we 
have 

n{Vi u V2) = n{Vi) ^n{V2) . (13) 

Now, let the density matrix of the system be p. We can define a density matrix pv for 
a volume V by tracing p over l-i{—V). Its entropy is SiV) = — Trpy log(py). Thus, the 
density matrices corresponding to the subsets satisfy the compatibility condition 
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Pv — Trn{v')Pvuv' , 



(14) 



where V is any set that has measure zero intersection with V. A continuous quantum system 
is defined here as a set of Hilbert spaces and density matrices corresponding to subsets of 
R*^ satisfying the conditions (13) and (14). The strong subadditive inequahties then take a 
geometric form given by (10) and (11), where now A and B represent subsets of R*^ [15]. 

The problem of building a continuous quantum system can be presented in a slightly 
different way where the complete system is constructed out of the finite volume subsystems. 
These may be defined for example by giving a Hamiltonian and boundary conditions, and 
taking a finite volume density matrix pv- Then, to glue all elements into a unique system 
we have to impose (13) and (14) on the local Hilbert spaces and states. 

When the total state p is invariant under some transformation U in R'^ the entropy 
S{U{A)) is equal to S{A). This symmetry combined with the properties of the entropy can 
strongly constrain the function S. We say that a positive function S on subsets of R'' is an 
Euclidean entropy function if it is invariant under rotations and translations, and satisfies 
SSA and SSB. 

Relativistic setting 

In the relativistic case a pair of time-like related subsets contains degrees of freedom 
that are not mutually independent. A scheme as the previous one for an Euclidean system, 
where (13) and (14) hold, must correspond to every Cauchy surface in spacetime. 

All previous work on geometric entropy assumes that we can divide a Cauchy surface 
for the spacetime into pieces and trace the total state over the degrees of freedom on one of 
them to obtain a density matrix. This assumption is also made here, since it is fundamental 
to the very definition of a geometric entropy. 

However, this is at the root of the divergences encountered in quantum field theory 
calculations. The complete localization of states in a given region is forbidden by standard 
axioms in QFT (see [17] for a clear discussion on this point). The covariant regularizations 
that deal with the infinites arising in perturbation theory can not in general avoid the 
divergence in the entropy, which shows up at the level of free fields. Also, this is not softened 
by supersymmetry. Thus, as these divergences suggest, there are two types of enigma related 
to the black hole entropy formula. The first is given by the proportionality between area 
and entropy, and the second one is that the entropy is surprisingly finite. Here we do not 
deal with this last and probably deeper problem. It is possible that a finite entropy would 
require quantum gravity input (see also a different idea in [18]). However, the existence of 
a good covariant definition for a finite geometric entropy in the context of QFT can not 
be discarded. But, in that case, the role of the set boundaries can not be taken as purely 
passive. In other words, to define a bounded system Hilbert space, boundary conditions are 
needed. The precise compatibility conditions that these bounded systems should have in 
order to be used for computing a geometric entropy are given by (13) and (14). This will 
be enough characterization for our purposes. Any non ambiguous definition for a localized 
entropy beyond the classical limit should suffer from the same problem. In particular, this 
is the case for the various entropy bounds present in the literature, which are suggested by 
ideas related to the physics of black holes [19]. In the semiclassical level they should take 
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(a) 



(b) 



FIG. 1. (a)- A causally closed set in 1 + 1 dimensions. It is the causal development of the 
spatial surfaces A and A'. We say that A and A' are Cauchy surfaces for A. Both of these sets 
can be continued to a Cauchy surface for the Minkowski space using the same spatial set Ci. The 
marked points on the left and right corners of the diamond shaped set A represent the spatial corner 
of A. (b)- Two commuting causally closed sets in 2 + 1 represented here by their Cauchy surfaces 
representatives A and B. These surfaces intersect in the dotted line in this example, dividing A 
into Ai and A2 and B into Bi and B2. The commutation imposes that A2 B and B2 C A. 
Note that the spatial corners of A and B, which are the boundaries of A and B respectively, are 
spatial to each other. The causally closed sets Ay B and A/\B that enter into the SSA relation 
are generated by A\ U Bi and A2 U B2 respectively. 

into account the geometric entropy of quantum origin, what can make them much more 
constraining and interesting [20]. 

Therefore, here we assign a density matrix and an entropy to any subset ^4 of a Cauchy 
surface C. In the following we restrict the analysis to the A with compact closure. The 
conditions (13) and (14) are assumed to hold among the subsets of any such surface C. 

The causal development A of A is the set where the solution of a wave like equation is 
uniquely determined by the initial data on A. More generally, for any set O not necessarily 
space-like, the causal development is the set of points x such that every inextensible time-like 
curve through x cuts O. We often say, abusing of the terminology, that a space-like set A is a 
Cauchy surface for A. Consider two subsets A and A' of different Cauchy surfaces C and C, 
but having the same causal development A = A' (see fig.(l)a). We can extend A and A' with 
the same set to form a global Cauchy surface, that is, we can choose Ci = C — A = C — A' . 
The unitarity of the causal evolution implies that we have 

S{A) = S{A') = S{A) . (15) 

This means that there is no loss nor gain of information in passing from A to A'. In the 
Heisenberg picture both p{A) and p{A') are obtained by tracing the total state over the 
same degrees of freedom in Ci = C — A = C — A'. 

Then, A can be seen as an equivalence class of pieces of Cauchy surfaces with the same 
causal development, density matrix and entropy. Note that all Cauchy surfaces for A share 
the same boundary, the spatial corner of A (the subset of the closure of A that is not time- 
like related to any point in A, see the figure (1)). The sets of the type A are called here 
causally closed, and they constitute the natural domain of S in the relativistic case. We 



7 



often make no distinction between the set A and an element A on the equivalence class it 
represents, and denote them with the same capital letter. 

Now we identify the conditions allowing to apply the strong subadditive property for any 
two given causally closed sets A and B. As a consequence of eqs. (13) and (14) the SSA and 
SSB inequahties in the form (10) and (11) take place for the entropies of two subsets A and 
B of the same global Cauchy surface C. In that case the sets An B and AU B appearing in 
the smaller side of the strong subadditive relation are Cauchy surfaces for AClB and the the 
causal development of AUB respectively. However, it is possible to apply SSA and SSB even 
in some cases where A and B do not belong to the same C. As the entropies really depend on 
the equivalence classes A and B, this happens if there are different Cauchy surfaces A' and 
B' for A and B that belong to the same global Cauchy surface C. We say that two causally 
closed sets A and B (or any pair of Cauchy surfaces for them) commute if there is at least a 
pair of the respective Cauchy surface representatives that belong to the same global Cauchy 
surface (see fig.(l)b). It can be seen that an equivalent condition for the commutativity of 
A and B is that their spatial corners are spatial to each other. In particular, if A includes B 
or they are spatially separated, they commute. The inequalities SSA or SSB can be applied 
only for the entropies of commuting elements. For these, the right hand side of (10) and 
(11) has to be calculated with any pair of Cauchy surface representatives belonging to the 
same global Cauchy surface. 

There is a more elegant way to define the concept of causally closed sets and commutativ- 
ity that has interesting connections with quantum mechanics [21]. Write a; ~ y for two points 
that arc time-like related. Given any subset O of the spacetime define its causal opposite 

as the subset of points which are spatial to O, O"*" = {x/x y,x ^ y, for every y G O}. 
The causally closed sets are just the subsets O which coincide with their double opposite, 
O^^ — O ^. Curiously, this sets form an orthomodular lattice (also called quantum logic). 
This structure is also shared by the physical propositions in quantum mechanics, or equiva- 
lently, the orthogonal projectors on the Hilbert space. In any orthomodular lattice there are 
three operations, the opposite, the meet A A B and the join A\/ B of two elements, and an 
order relation C. In the present context they are given respectively by the above opposite, 
the intersection Ad B, the double opposite of the union {A \J B)-^-^, and the set inclusion 
relation C. There is a definition of commutativity of two elements for any orthomodular 
lattice that coincides with the usual commutativity of operators when applied to projectors 
in the Hilbert space. Two elements A and B commute when 

A^(AAB)y (AaB^). (16) 

This definition coincides with the one given above for the commutativity of causally closed 
sets. Thus, the compatibility requirement that two causally closed sets have Cauchy surfaces 
which can be extended to a common global one can be put in the same mathematical terms 
as the compatibility condition for two projection operators to be simultaneously measurable. 
With this notation, the SSA and SSB relations for two commuting sets A and B are 



This equivalence holds for bounded sets. The causally closed sets O = O are exactly the 
domain of dependence of achronal sets bounded in time [22] . 
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S{A) + S{B) > S{A VB) + S{A A B) , 
S{A) + S{B) > S{A A B^) + S{B A A^) , 



(17) 
(18) 



respectively. 

A relativistic entropy function is defined as a positive^ function on the causally closed 
sets in Minkowski space, symmetric under Poincare transformations, and satisfying SSA and 
SSB for commuting pairs of subsets. The use of the relation SSB in the proofs below may 
be avoided, but it considerably simplify the demonstrations. 

Thus, the problem of finding the functional dependence of the geometric entropy takes 
a precise mathematical form, being equivalent to find the most general relativistic entropy 
function. We note that a linear combination of solutions with positive coefficients is also a 
solution. 

A similar problem can be stated in any space-time, but in absence of the Poincare 
symmetry the number of solutions can be very big. For example, the flux of any future 
directed conserved current through a Cauchy surface for a causally closed set A is always 
a solution of SSA and SSB since it satisfies the strong subadditive equation (rather than 
inequality). Of course this solution is not Poincare invariant in Minkowski space. Symmetries 
with orbits that are exclusively time-like (or space-like) do not impede this type of extensive 
solutions. 

In this paper we are only using causally closed sets with compact closure, and, as in 
the Euclidean case, make no distinction between sets differing by zero measure subsets. 
The game to compare the entropies on different sets using strong subadditivity is more 
subtle than in an ordinary additive measure theory. Because of that we further restrict 
the class of subsets. To see what are the difficulties about, consider first the polyhedra on 
the Euchdean space R*^. The intersections and unions of polyhedra are again polyhedra. 
Thus, the relations SSA and SSB can not be used to obtain information for the entropy 
of sets with curved borders from the entropy on polyhedra. In order to focus on the most 
restrictive relations coming from the entropy properties and the symmetries, in the following 
we consider the domain of the Euclidean entropy functions to be the class of all polyhedra. 
Note, however, that these can have an arbitrarily large number of faces as small as we 
want, and, in this sense, we can approach any set by polyhedra. To connect the results 
about polyhedra with other kind of sets some continuity condition is needed, and we do 
not deal with this problem here. In the relativistic d + 1 dimensional case we consider the 
class of relativistic polyhedra, that is, spatial d dimensional sets (or their corresponding 
d + 1 dimensional causally closed sets) formed by union of a finite number of polyhedra on 
spatial hyperplanes. These also form a closed class under union and intersection when taken 
between commuting sets. 

A note on the terminology. We call vol(X) the volume of a dimensional Euclidean 
polyhedron X, and area(X) to the area of the d — 1 dimensional boundary of X. For a 
relativistic d + 1 dimensional polyhedron X we call area(X) the area of its spatial corner, 
which is the boundary of any Cauchy surface for X. A relativistic entropy function S 



^There was an indication of a negative geometric entropy when gauge fields are present [6]. This 
troubling result was finally corrected [8,10]. 
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on d + 1 dimensions when restricted to polyhedra on a unique hyperplane gives place to 
an Euclidean entropy function on d dimensions. This latter is independent of the chosen 
hyperplane. When it is possible, we identify these functions, and call them with the same 
letter. 

IV. ENTROPY IN AN EUCLIDEAN SYMMETRIC SYSTEM 
A. The one dimensional case 

Consider a system on the real line and a translational invariant density matrix p. The 
connected compact sets on are intervals, determined up to a translation by their length 
X. On these intervals the entropy can be written as a function on a real variable S{x). The 
results we obtain here about S{x) can be found for example in the ref. [15]. 

Take two intervals of length {x + y)/2, with x and y positive and y > x, and arrange them 
so that their intersection has length x and consequently their union has length y. Applying 
strong subadditivity we have 

S (^) > ls{x) + ls(y) . (19) 

Suppose now that S{x) is infinite. Then it follows from (19) that S{z) must be infinite for 
{x/2) < z < oo. Repeatedly using this argument we have that the function S{x) is either 
finite for all x or always infinite. 

Similar arguments can be used to show that all the entropy functions considered in the 
paper are either infinite or always finite. We are not coming back to this point and assume 
a finite S in the following. 

The relation (19) for the function S{x) is called weak concavity, while the concavity is 
the property 

S{Xx + (1 - X)y) > XS{x) + (1 - X)S{y) , (20) 

for any < A < 1. A function is concave when the segment joining any two points in its graph 
lies below the function curve (see the fig. (2)). Weak concavity is in fact very near concavity, 
since applying (19) repeatedly it is easy to obtain eq.(20) for any rational A with a power of 
2 as a denominator. For example, we have S{3/4:X + l/4i/) > jS{x) + jS{l/2x + l/2y) > 
^S{x) + \S{y). To show that S{x) is actually concave suppose that x < z < y and, violating 
concavity, 

Siz) = ^S(x) + '-^S{y) - 5 , (21) 
y — X y — X 

with 5 > 0. Let w = Xx + [1 — X)y with A e [0, 1] a rational with a power of 2 as 
a denominator. We can take w = z + e with e a positive number as small as we want. 
Positioning two intervals of length z such that their union has length z + e = w, and their 
intersection has length 2; — e, we have from SSA and (21) 

S{z - e) < S{z) - 5 + 0{e) . (22) 
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X 

FIG. 2. A non decreasing concave function S(x). Two points on the curve S{x) determine a 
segment that hes below the function graph. The hmit of S{x) for x ^ is 7, and the curve slope 
approaches to s as x goes to infinity (however, S{x) does not necessarily approach any straight line 
asymptotically) . 



Thus, as near to z as we want there are points in which, let say, S < S{z) — |5. By repeating 
this argument for a sufficient number of times we could arrive at the result that there would 
be points as near to z as we want that have negative S. As the entropy is positive, equation 
(21) must be wrong, and the function S is concave, equation (20) holds for any A e [0, 1]. 
The converse is also true, concavity implies strong subadditivity for intervals. 

Suppose that for some y > x it is S{y) < S{x). Then the function S is decreasing for 
any z > y, because a consequence of concavity is that the slope of the function S always 
decreases. This would lead to S{z) < for high enough z. To avoid it the function S 
must be non decreasing everywhere. The monotonicity of the entropy appears here as a 
consequence of translation symmetry. However, order by inclusion is not warranted for the 
entropies of sets formed by several disjoint intervals, or connected sets in more dimensions. 
The concavity also imply that S{x) is continuous. 

Being nondecreasing and positive S{x) has a well defined limit as x goes to zero. We 
call 7 this limit. Given x < y we have from concavity that 

S{x)>-f + x^i^^^. (23) 

y 

Since 7 > it follows that S{x)/x isa, non increasing function. This leads to the existence of 
the limit limj-^oo S{x) jx = s, with s > 0. Therefore there exist a well defined notion of mean 
entropy in the system (otherwise S{x)/x could increase to infinity or oscillate indefinitely). 
For big enough sets the entropy is approximately extensive, at least if s 7^ 0. 



B. The limit for the entropy of small one dimensional sets 



The measure zero sets have zero entropy. However, sequences of sets with measures 
converging to zero can have non zero entropy limit. The limit of S{x) when x goes to 
zero is in general a positive number 7. Consider now the case of a set formed by two 
intervals of lengths xi and X2 separated by a distance yi. We are interested in the limit 
Z = Xi + ?/i + X2 — >^ 0. If e is a small positive number we can choose / such that S{1) < 7 + e. 
Placing an interval of length y[ <l,y[> yi as shown in fig.(3)a it follows from SSA that 
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1 I 1 

yi X2 ^1 yi X2 

(a) (b) 

Xn+l ''^n+ Xi 



^1 Xn Xj Xn Xn+l 

(c) (d) 

FIG. 3. Geometric constructions used to prove the formula (28) for the entropies of one dimen- 
sional sets in the limit of small size. 

S{xi, y,, X2) + S{y[) > S{x[,yi, x'^) + S{1) , (24) 

with x[ < xi, X2 < X2, x[ + X2 + yi = y'l- Taking into account the monotonicity of S{x) we 
have 

S{xi, yi, X2) > S{x[, yi, x'2) , (25) 

for any x'l < xi and x'2 < X2- Similarly, from the same construction of fig. (3) a but now using 
SSB we have 

S{xi,yi,X2) > S{x[,y[,X2) - e, (26) 

for any y[> yi, x[ < xi, X2 < X2, and Xi+yi+X2 = x[ +y[ + X2. Thus, reducing the length 
of the two intervals reduces the entropy (neglecting contributions of order e). The converse 
is given by the construction of figure (3)b. Applying strong subadditivity we get 

S{xi,yi,X2) > S{x[,y[,X2) - e, (27) 

for any x[ > Xi, y[ < yi, and Xi + yi = x[ + y[. 

Prom (25), (26) and (27) we see that S{xi, yi, X2) converge to a limit 5 when xi +yi + X2 
go to zero, and this limit is independent of the way in which the interval lengths Xi and X2 
go to zero. 

Now consider the entropy S{xi,yi, ...,ym^i,Xm) of a set formed by m disjoint intervals 
of lengths Xi, ...,Xm separated by distances yi, ^m-i- We analyze the entropy in the limit 
when I — xi + yi + ... + Xm go to zero and prove by induction that this depends only on the 
number of connected components m, and is given by 



where 



\imS{xi,yi,...,yrn-i,Xrn) = Sm = a + pm, (28) 



a = 27 -5, (29) 
/? = 5-7. (30) 
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For m = 1 and m = 2 we have proved this formula. Suppose it is vahd for m < n. Prom 
the constructions of fig.(3)c and fig.(3)d and SSA we have 

S{xi,yi, Un-l, Xn) + S{Xn, y„, Xn+l) > S{xi, t/i, y„, Xn+l) + S{Xn) , (31) 
S{xi, yi, y„, Xn+l) +S{Xn + yn + Xn+l) > (32) 

S{Xn, y„, Xn+l) + S{xi, yi, yn-1, Xn + yn + Xn+l) , (33) 

respectively. Taking the limit of small xi + |/i + . . . + y„ + Xn+i in these inequalities and using 
the induction hypothesis it follows that (28) is also valid for m = n + 1. Then it holds for 
any integer m > 1. 

It is immediate from subadditivity and from (28) that 

< 7 < 5 < 27 . (34) 

These relations for the parameters simply translates into the positivity of the independent 
variables a and /3, which do not have to satisfy any constrain. We will also write (28) as 

Sn,^^ + {m-l)(3, (35) 

where 7 > /3 > 0. 



C. The multidimensional case 



Consider now a translation invariant state in R''. Let the unit vectors Vi, ...,Vd form an 
orthogonal basis of R'^. Define the family of rectangular polyhedra as all the translates of 
the R{ai, .... a^) = {x/x = Yli=i ^iVi, Aj G [0, Oj]}, where (ai, .... a^) are the sides lengths. We 
often use just a symbol R for these polyhedra, what should not be confused with the symbol 
for the reals R. The entropy can be written as a function on d variables S{ai, ...,ad) on 
these sets. The same constructions made above for the one dimensional case apply here on 
each direction separately. Therefore S{ai, ...,aa) is a non decreasing and concave function 
in each a^. Then the entropy is ordered by inclusion among the rectangular polyhedra. 
Also, by concavity S{{ai, ad)/{ai x ... x aa) is non increasing in each separately. Then 
the entropy divided the d dimensional volume for has a limit for the rectangular polyhedra 
when the volume goes to infinity along any given increasing sequence When all the 
sides of Rn go to infinity, the limit s of entropy over volume must be the same number 
independently of the particular sequence. This is because giving two of such sequences of 
rectangular polyhedra with sides going to infinity, every element of one sequence must be 
smaller than some element of the other sequence, and in consequence none of the limits of 
the entropy over volume could be bigger than the other. Thus s is really the entropy per 
unit d dimensional volume of the system. If the volume goes to infinity but one of the sides 
remains bounded the limit of entropy per unit volume could differ from s. Physically, this 
means that the effects of the boundary may always be relevant in this case. 

When in addition to symmetry under translations there is rotational symmetry this result 
can be generalized to different sequences of polyhedra. Any sequence where the volume goes 
to infinity at a greater rate than the bounding area (going to infinity in the sense of Van 
Hove), and where the number of faces remains bounded, has the limit s of entropy over 
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FIG. 4. Eating from inside a strip D to a smaller strip D' (the right side of D' is the dashed 
line), by means oi A — B. One has to position A inside D (shown with the dotted line) and take out 
from D the sector corresponding to A — B = C. Then translating down A inside D and repeating 
the operation it is possible to eat all D — D'. 

volume, which is independent of the particular sequence [15]. It is in this type of sequences 
that one expects the surface terms contribution to the entropy to be negligible in the limit. 
We do not use this result here. Instead we are more interested in the properties of the 
entropy function for sets with finite area and vanishing volume than in the infinite volume 
limit. These flat sets may not be so interesting in statistical mechanics but they are crucial 
for the relativistic geometric entropy. The boost symmetry allows to amplify their width 
and reduce most of the general relativistic problem to the case of flat Euclidean sets. 



D. Entropy for flat sets 

Let us consider polyhedra on R'^ of the form X x [xi,yi, ...,ym-i, Xm], where X is a d — 1 
dimensional polyedron and [xi,yi, ym-i, Xm] is a set with m connected components in R^. 
Here Xi, ...,Xm are the single component lengths and yi, ym-i are the separation distances 
between adjacent components. Call the corresponding entropies S{X,xi,yi, ...,Xm)- For a 
fixed X these polyhedra form a translation invariant one dimensional system, and we have 
from equation (35) that 

hm SiX,xi,yi,...,Xrn)^liX) + im-l)f3{X), (36) 

Xl+yi + ...+Xm^O 

where '~f{X) and P{X) are positive and Euclidean invariant functions on polyhedra in R'^"^, 
and 7(X) > f3{X). It follows from SSA and SSB in d dimensions that 7(X) and l3{X) also 
satisfy SSA and SSB. Note that we use the same symbols /5 and 7 that represent constants 
in the one dimensional case for denoting functions in the multidimensional case. 

In addition, these Euclidean entropy functions are also ordered by inclusion. To show it 
we first introduce a method that we can call the 'eating from inside' procedure. This was 
developed in [13] as a tool to prove certain monotonicity relations. 

Suppose we have two sets A and B with B C A, and we want to obtain information about 
the difference of their entropies. Call C — A — B. At our convenience we choose a set D with 
Ac D. Then from strong subadditivity it follows that S{A) + S{D - C) > S{D) + S{B). 
This can be written as S{A) — S{B) > S{D) — S{D — C). Now, if a symmetric copy Ai of 
A can be situated inside D — C we can repeat the procedure with D — C in the place of C, 
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FIG. 5. The geometric construction used to show that the function (3 is proportional to the 
volume. The solid line represents a transversal view of the two component set A. Drawn with 
dashed lines are the two component set Bi on the left and the two component set B2 on the right. 
All the connected components of these sets have the same height a. 

leading to S{A) - S{B) > S{D - C) - S{D -C-Ci), where Ci is the symmetric copy of C 
corresponding to the same transformation that carries A into Ai. Summing up this and the 
previous inequality we have that 2{S{A) — S{B)) > S{D) — S{D — C — Ci) . The procedure 
can be continued if at the n*^ step we can place a copy of A inside what is left from D after 
subtracting all the previous symmetric the copies of C (see fig. (4)). In that case we have 

S{A)-S{B)>-{S{D)-S{D')), (37) 
n 

where D' = D — C — Ci — — C„_i. In some situations the right-hand side in (37) is known, 
specially if the limit n — > 00 can be taken. Prom this we can obtain information about the 
difference of entropies of two sets A and B ordered by inclusion. We say that to obtain (37) 
we cat from inside D to D' by means oi A — B. 

Take now any d — 1 dimensional polyhedra X and Y with V C X, and let the one 
dimensional set Z = [xi, yi, ...,ym-i, ^m] = [ct; 0,, cl-i o\ be formed by m intervals of length 
a, with separation distance a between adjacent intervals. Consider the sets A — X x Z, 
B^YxZ,C^{X-Y)xZ, and D ^ X x [2maN] . Here [I] means a single interval of 
length /, and N is a. big integer. Then it is easy to see that translating A, B, and C on the 
direction perpendicular to X and Y we can eat from inside 2N copies of C from D, leading 
to the remaining set D' = Y x [2maN]. Application of (37) gives 

777/7 

The results for one dimensional translation invariant systems imply that S(U x [I]) /I has a 
limit su for / going to infinity. This implies taking N ^ 00 that 

S(A) - S(B) > ma (sx - sy) ■ (39) 

Now, taking the limit a — * we have from this inequality that '~f{X) + (m — l)l3{X) is 
ordered by inclusion for any m > 1. Consequently 7(X) and P{X) are ordered by inclusion. 

Moreover, the functional form of P{X) can be calculated explicitly. Let us take A — 
X X Z with Z = [xi.yi, X2] = [a, a/2, a] and X a d — 1 dimensional polyhedron. Also take 
Bi — Yi X Z and B2 — Y2 x Z with Yi and Y2 two d — 1 dimensional polyhedra such that 
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X = Yi U Y2 and Fi fl I2 = 0- Suitably placing Bi and B2 to have empty (or measure 
zero) intersection and to cover the gaps on A (see fig. (5)), we can arrange that the union 
AU BiU B2 = X X [(13/4)a] and the intersections A n Bi = Yi x Z' , A n B2 = Y2 x Z', 
where Z' — [xi,yi,X2,y2,X3] = [a/4, a/2, a/4, a/2, a/4]. Then, from strong subadditivity for 
the pairs {A,Bi) and {A\J Bi, B2), and taking the hmit a ^ we get 

(3{X) > (3{Y,) + (3{Y2) . (40) 

This inequality is opposite to the subadditive inequahty for /3 since X — Yi[JY2. Thus the 
equation holds, 

P{X) = P{Y,) + ^{Y2) . (41) 

With A^'^^^ copies of a rectangular polyhedra R with sides (/i, Id-i) we can exactly cover 
a rectangular polyhedron U with sides {N li, ...,Nld-i). Using (41) it is 

The fraction on the right hand side converges to a limit 2s for large N. The factor 2 is 
chosen for future convenience. Thus, we have that the function /3 is proportional to the 

volume for any rectangular polyhedron. To show that f3{X) is proportional to the volume 
for any polyhedron X divide the space R^^^ with a cubic grid formed by copies of the 
cubic element A. Call r~{A,X) the union of elements of the grid that are included in X, 
and r~^{A,X) the union of elements of the grid having non empty intersection with X. 
Using (41) and that /3 is proportional to the volume on cubes it follows that the entropies 
S{r^{A, X)) = 2s vol(r^(A, X)) are proportional to the volumes. Then, using the order 
by inclusion of (3, the fact that yo\{T-{A,X)) < vo\{X) < vol(r+(^, X)), and that these 
volumes coincide in the limit of small A, we obtain 

(3{X) ^2svo[(X) (43) 

for every X. Note that the only symmetry we used was translational invariance. 

The function a = 7 — /3 is positive and strong subadditive because P{X) — 2svol{X) 
satisfies the strong additive equality, P{A) + (3{B) = (5{A n S) + (5{A U B). 

Summarizing, we have shown the following theorem on the entropy limit for fiat sets 

Theorem 1: For an Euclidean invariant d dimensional system the limit of the entropy 
for sets of the type X x [xi, yi, x^] as all the xi, yi, go to zero exists and is given 

by 

lim 5(X,xi,yi,...,xJ=7(X) + (m-l)2svol(X) , (44) 

where 7(X) is an ordered by inclusion d — 1 dimensional Euclidean entropy function, and 
7(X) > 2svol(X). 

One could wonder about the entropy limits for lower dimensional objects. Define 7 and 
(3 as the operators that transform Euclidean d dimensional entropy functions into Euclidean 
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FIG. 6. Two commuting sets with straight Cauchy surface size V xx' in 1+1 dimensions intersect 
in another diamond-hke set of size x', and their union is, after causal completion, a diamond set 
of base size x. 

d — 1 dimensional entropy functions and whose action is defined as 75 = 7(5*) = 7, /35 = 
I3{S) = f3, with 7 and f3 as above. The monotonicity of 7 and f3 imply that /3 o 7 = 0, and 
/3 o /3 = 0, and the explicit form of /5 leads to 7 o /5 = 0. However, the composition 707 
can be different from zero when apphed to a general entropy function S. We call 75 for 
n — 0, d the composition of S with n times the operator 7, and define 75 = S. 

V. GEOMETRIC ENTROPY IN MINKOWSKI SPACE 
A. The 1 + 1 dimensional case 

Before embarking in the calculation of the functional form of the entropy in the general 
d + 1 dimensional Minkowski space we focus on the simpler 1 + 1 case. This allows to 
sec without the complications of higher dimensions the important role played by the boost 
symmetries. The result for the entropy and the method to obtain it are similar to the ones 
for the 1 dimensional Euclidean small sets. 

The connected causally closed sets are determined in 1 + 1 dimensions by a single number 
X, the size of the unique straight Cauchy surface (the size of the diamond base, see fig. (6)). 
The corresponding entropy S{x) is a positive, non decreasing and concave function because 
of translation symmetry. As above, call 7 the constant \imx-ioS{x). 

The additional boost symmetry greatly constrain the function S{x). To see this, consider 
the construction of fig. (6). Given x' < x we place two commuting boosted one component 
sets of size y/xx' in such a way their intersection is a one component set of size x' and their 
union is another one component set of size x. The application of SSA gives 

S{x) + S{x') < 2 S{V^) . (45) 

Making x' as small as we want in this inequality we prove that S{x) < 7. As S{x) is non 
decreasing this means that S{x) = 7 is constant. Due to the boost symmetry we can form 
a Cauchy surface for a given diamond D using two boosted diamonds of size as small as we 
want disposed along the null surfaces of D. This, and the subadditive law, is what impedes 
an entropy increasing like the volume. 
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(a) 
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0: 
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(c) 



(d) 



FIG. 7. Two dimensional constructions used to prove Theorem 2. In (a), (c) and (d) it is 
considered the SSA between the sohd hne and the dashed hne set. In (b) both the SSA and the 
SSB between the solid hne set and the central dashed diamond are used. 

The entropy for sets consisting on several components spatially separated can in principle 
depend on the lengths of the individual diamond bases, their separation distances, and 
relative angles. 

Let us start with the case of two component sets. In the construction of fig.(7)a the union 
of a two component set A with a single diamond S is a set A', having two components, one 
of which is bigger than the corresponding component in A. Application of SSA gives 



Thus, the entropy decreases with increasing size of the individual connected components. 
Now, the application of SSA and SSB to the construction of fig.(7)b leads to the inequalities 



where A' and A" are both included in A but A' shares with A the internal spatial corners 
and A" shares with A the external spatial corners. Combining these relations we can prove 
that the entropy of any two component set X is greater than the entropy of any other two 
component set X' where each individual component of X' is included in a corresponding 
component of X. Thus, the entropy of the two component sets is increasing and decreasing 
with inclusion, and then it must be a constant. Call 6 this constant. 

For a greater number of connected components the entropy can be expressed in terms of 
7 and 5 and we show the following 

Theorem 2: The most general form for a relativistic entropy function in 1 + 1 dimen- 
sional Minkowski space is given by 



S{A)+'y>S{A')+'y. 



(46) 



S{A) + 7 > S{A') + 7 , 
SiA)+^>SiA")+j, 



(47) 
(48) 



S{Z) = Sm^l+{m-l)(3, 



(49) 
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where m is the number of connected components of Z and 7 > /3 > 0. 

In terms of 5 it is /3 = 5 — 7. The proof is very similar to the corresponding one for 
the small one dimensional sets in Section 111. To show (49) by induction on the number of 
components note that for m = 1 and m = 2 this equation holds. Supposing it is valid for 
m < n we have from the constructions of fig.(7)c and fig.(7)d 

Sn+liX)+^>Sn + 5, (50) 
Sn + S>Sn+l{X)+l. (51) 

This proves the formula (49). It follows from positivity and subadditivity that the parameters 
have to satisfy the relations 

0<7<5<27, (52) 
< /? < 7 . (53) 

The meaning of these constraints is clarified if we write as a sum 

Oi + pm, (54) 

of two independent solutions with independent and arbitrary positive coefficients. Here 
a — '-f — 0. One of the solutions is the constant term and the other is topological, being 
proportional to the number of connected components m. It is easy to check that (54) does 
indeed satisfy SSA and SSB for two arbitrary commuting sets. Thus, this formula, or eq. 
(49), represents the most general solution for the entropy in 1 + 1 dimensional Minkowski 
space. 



B. Multidimensional case d> 1 

In this Section we show the theorem 

Theorem 3: The most general form for a relativistic entropy function on the relativistic 
polyhedra on d + 1 dimensional Minkowski space, with d > 1, is given by 

S{X) = ao + s area(X) , (55) 

where ckq and s are non negative constants and area(X) is the area of the spatial corner of 
X. 

We divide the proof in several parts. 



1. Relativistic Lemmas 

To apply SSA or SSB between two sets we have to ensure that they commute, what 
complicates the proofs for d > 1. In the following we incorporate into the Euchdean function 
given by restricting the relativistic entropy function to a single spatial hyperplane some 
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information coming from the boost symmetry. With these additional properties at hand we 
will be able to work in Euclidean space where all sets commute. 

An example where commutativity fails is given by adding one perpendicular spatial 
dimension to the sets in figure (6). The two boosted intervals of length \/xx' there can be 
seen as the transversal section of two boosted rectangles in 2 + 1 dimensional Minkowski 
space. However, these rectangles do not commute when they are in such a position since 
their boundaries are not spatial to each other. 

There is a limit where commutativity can be applied in an effective form. Consider 
the product x Z of a big d — \ dimensional rectangular polyhedron R and a relativistic 
polyhedron Z in 1 + 1. We mean here by the product of a polyhedron R with a relativistic 
polyhedron Z the relativistic polyhedron generated by the Cartesian product of R with a 
Cauchy surface for Z, since the Cartesian product of R with the relativistic polyhedron 
corresponding to Z is not a relativistic polyhedron on o? + 1. Given two commuting 1 + 1 
dimensional sets Z and Z' the sets RxZ and RxZ' may not commute. However, taking two 
fixed commuting Z and Z', we can construct a relativistic polyhedron U commuting with 
Rx Z hy deforming RxZ' near the borders of R. More precisely, we take U as the domain 
of dependence of the union of {R x Z') r\{RxZ) and (i? x Z') n{Rx Z)^ . Here is the set 
formed by all the points spatially separated from A. When the sides of R are big, U differs 
from i? X by a set located along the boundary of i?, and whose transversal size remains 
fixed if the sizes of R are increased. Thus, the difference in entropy between RxZ' and U 
can be bounded by a number proportional to the border area of R by using intersections 
and unions of rectangular polyhedra with bounded sides. If we consider the limit of the 
entropy S{R x Z) divided by the volume vol(i?) for the sets of the type RxZ with R going 
to infinity the terms growing like the border area of R can be neglected. Therefore, in this 
limit we can use commutativity between the sets RxZ and RxZ' with commuting Z and 
Z', and use the two dimensional results. We have the 

Lemma 1: Let it! be a d — 1 dimensional rectangular polyhedron and Z a 1 + 1 relativistic 
polyhedron then 

S(Rx Z) , , 

i™ Lum = Too + 2s m - 1 , 56 

where m is the number of connected components in Z, s is defined by Theorem 1, and 7oo, 
with > 2s, is the limit 7(i?)/vol(-R) for R going to infinity. 

The explicit form of the right hand side follows from Theorem 1, since the two dimensional 
solution tell us that this term must be a constant plus a constant times the number of 
components of Z. Then we can cahbrate the parameters from the special case of Z lying in 
a single spatial line, Z — [x, y, y, x], for small x, y. 

In particular the Lemma applies to the sets Z = [xi,yi, ...,Xm] that lie in a single spatial 
plane, with arbitrary Xi,yi, x^- Thus, it removes the constrain of small Xi+yi + ... + x^ 
in Theorem 1, in the limit of R going to infinity. Of course, this includes information about 
boost symmetry and it is not the case for a general Euchdean function. 

With the notation of Section IV consider now the entropies S{R,Xi,y,X2) for the sets 
X = Rx [xi,y,X2\ contained in a single spatial hyperplane, and where i? is a fixed d — 1 
dimensional rectangular polyhedron. Call Xi — R x [xi] and X2 — R x [2:2]. Theorem 
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FIG. 8. Two dimensional cut of the geometrical construction used to prove inequality (58). 
The lines at 45° are null surfaces. The dotted lines represent the spatial hyperplanes containing 
the defining Cauchy surfaces for Y = Y1UY2 and W = Wi U W2- 

1 showed the form of this function for the hmit in which all xi,y,X2 go to zero. We can 
do better in the relativistic case, and find the limit for xi and X2 going to zero and any 
y. The first inequality we obtain comes eating from inside U x [xi,y,X2], where f/ is a big 
rectangular polyhedron, to U x [xi] by means oi X — Xi. Using the Lemma 1 and Theorem 
1 this leads to 

lim S{R X [xi, y, 2:2] ) > -f{R) + 2s vol(i?) . (57) 

Then consider the following construction shown in the figure (8) through a two dimensional 
projection. The set y = Yi U 1^2 is of the form R' x [x\,y',x'2\, but it appears highly 
boosted with respect to X in the direction perpendicular to R. The same happens for 
W — W1UW2 = R" X [x'l, y", X2]. The rectangular polyhedron R" Q R' Q R, the transversal 

), and the boost factors are chosen in such a way that Yi C Xi, 
Wi C Y2, W2 C X2, and Y2 is spatial to Xi and X2. We are taking then the limit of small 
(xi, X2, x[,y', x'2, x'{, y", X2) but keeping y fixed. In this hmit the sides of R' and R" can be 
taken to converge to those of R. Applying successively SSA to the pairs of sets {Xi,Y), 
(Xl \JY2,W), and (Xl \JY2UW2, X2) we have 

S{X U Y2) < 27(i?) - 7(i?") + 2s (vol(i?') + vol(i?")) + e . (58) 

for e as small as we want. Now, we eat from inside a three component set formed by the 
product of a big rectangular d — 1 dimensional set U and the relativistic 1 + 1 set given by 
the transversal projection of X U 12 (the set shown in figure (8)) to U x [xi, y, X2] by means 
of (X U Y2) — X. Taking into account Lemma 1 this leads to 

S{X U Y2) - S{X) > 2s vol{R') . (59) 

Combining the last two equations, taking a limit in which R", R', and R all have the same 
volume and shape, and using the continuity of 7 as a function of the sides of the rectangular 
polyhedron, we obtain the opposite inequality to (57). Then, we have 

lim S(Rx\xi,y,X2])^^(R) + 2svo\(R). (60) 
This equation allows us to prove the following 
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a 

FIG. 9. The entropy of X U Y2 is equal to S{X) + 2svol(i?) in the hmit of smaU a. The 
dotted Unes represent a two component strip-Uke set that can be eaten from inside down to a one 
component strip by means of {(X U Y2) — X). 

Lemma 2: Let X he a d dimensional polyhedron, R a, d — 1 dimensional rectangular 
polyhedron, Z — [a,y,a], Y = R x Z, and Yi and Y2 the connected components of Y, 

symmetric copies of R x [a]. Suppose that Y f] X = Yi and that there is an open set O 
including Y2 such that any translation of O in a plane parallel to R has no intersection with 
X (see fig. (9)). Then 

lim S{X U Y2) = S{X) + 2s wo\{R) . (61) 

a— >0 

The proof uses the previous result (60) and Lemma 1. Applying SSA to the pair [X, Y) 
we have 

lim Six U Y2) < S(X) + 2s vo\(R) . (62) 

a— >0 

The opposite inequality follows from eating from inside a two component large rectangular 
strip-like set, one component containing X, the other having transversal size equal to a and 
containing Y2, as show in fig. (9). 

2. Subtraction of the constant term 

Here we show that the limit of the entropy for rectangular polyhedra with vanishing 
small sides can be identified with the constant solution component of the general solution. 
In other words we show the following 

Proposition 1: Let S" be a relativistic entropy function. Then S — ao + S^, where 
is the limit value of the entropy for rectangular polyhedra of vanishing small sides, and 
is a relativistic entropy function with zero entropy limit for the rectangular polyhedra with 
vanishing small sides. 

We know that the general solution must have a constant term. However, the identification 
of it with cto and the subtraction of cto are not automatic, since we have to show that is 
still positive and satisfies SSA and SSB. An example where this is not the case is given by 
the two dimensional solution (49). If one subtracts the entropy limit for small intervals (7 
in the two dimensional case), the result is a positive but not strongly subadditive function. 
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(a) 



(b) 



FIG. 10. (a)- The entropies of two non intersecting sets X and Y satisfy the relation (65). This 

is shown by constructing X' = X U R and using SSA on the pair {X' , Y). (b)- For Y C X wc can 
show that S(X) + SiY) > S{X -Y)+ ao/2 by using SSB on the pair {X', Y), where X' = X U R 
and R is a smah rectangular polyhedron spatially separated from X. 

This is because subtracting 7 is different from subtracting the constant solution. There is a 
non constant solution, the one proportional to /3, that is constant over the one component 
sets and adds to 7. However, as it will be clear from what follows, in more dimensions 
there is no such solution. In particular, there are no positive topological solutions to the 
strong subadditive inequalities. Curiously, the Euler number E{X), which is the second 
solution appearing in the two dimensional formula (54), satisfies the strong additive equality, 
E{A) + E{B) = E{A U B) + E{A n B) in any dimension. However, for d > 1 it attains 
negative values. 

Now consider any relativistic polyhedron X. It contains a sufficiently small rectangular 
polyhedron R' lying in a spatial hyperplane V. We can generate a copy R of R' by translating 
R' on in a direction perpendicular to a phase F of R', and such that Rf] X — 0. For 
R sufficiently far from X there is a time-like hyperplane parallel to F in between X and R 
that divides the Minkowski space in two parts, one of which contains X and the other R. 
Let us call X' = X Li R (see fig. (10) a). Now, from Lemma 2 we have that for small enough 
R 

S{X') - e < S{X) < S{X') + e . (63) 

Then, by means of a tt angle rotation and a small translation on V we can make a copy of 
X' whose intersection with X' is a small rectangular polyhedron inside R. It follows from 
strong subadditivity that S{X') > (1/2) ao- Then taking the limit of small R we have 

S{X) > (64) 

for any X. 

Thus, we can safely subtract (1/2)q;o from the entropy without violating positivity. 
Strong subadditivity for the subtracted function will be automatic for the case where the two 
commuting sets X and Y have non empty intersection. When the intersection X (lY — $ 
we have to prove 

S{X) + S{Y) >S{XUY) + ^, (65) 

to show that the subtracted function is still subadditive. To see this, consider a set X' — 
X\JR s& above, but now we choose R included in Y (see fig(lO)a). Then 
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S{X') + S{Y) > S{X u y) + S{R) . 



(66) 



From (63) and taking limits the relation (65) follows. 

This shows that subtracting the constant solution ao/2 leads to a positive strongly sub- 
additive solution. It also satisfies SSB as follows from similar arguments based on the 
construction of fig(10)b for the case in which the inequality is applied to the sets X and Y 
with Y C X (otherwise there are two non empty sets on each side of the SSB inequality, 
which holds trivially for the subtracted function). 

Once we have subtracted q;o/2 we can repeat the procedure with the subtracted function 
S — ao/2 to subtract ag/A, and then ao/8, etc.. It follows that wc can subtract a constant 
solution xaQ for any < x < 1 and still get a solution. Thus S — is also positive and 
satisfies SSA and SSB. 



Taking into account the Proposition 1 we work with the subtracted function and assume 
CKo = 0. We show that ^{X) — /3{X) — 2svol(X), or, what is the same, a{X) — 0. 
To do so we need to make some geometrical constructions and use that the entropy for 
sets that are flat in more than one perpendicular direction can be safely neglected. In the 
notation of Section III this writes j"'{S) = for n = 2,...,d. It is possible to write this 
last equation and the result we are looking for, 7(X) = P{X) in a compact manner, that 
is, = (7 - /9) o 7g~^ = 0, for n = 1, d. Here we used that /3 o 7^"^ = for n > 2. 
Recall that 75 = S. Each is a function on the polyhedra on R''"", with ag — a and 
ag = ao = 0. 

Then, we show by induction that all the functions ag for n = l,...,d are zero in the 
relativistic case when ao = 0. For n — d this is true. We show that if a'p'^ — then = 0. 

Consider first a polyhedra T in R<^-"+i with a wedge-like form, T — {{xi, Xd-n+i)/^ < 
xi < 11X2+1', < X2 < 02, < Xd-n+i < cid-n+i}- Call R the d—n dimensional rectangular 
polyhedron with sides (02, ad-n+i)- We are interested in the limit of 75^^(T) as /i and u 
go to zero. Apply SSA and SSB to the constructions of fig(ll)a and fig(ll)b respectively, 
and using the induction hypothesis neglect sets that are small in more than one direction. 
We obtain two inequalities that are opposite in the limit of zero /i, u. They lead to the 
equation 



Thus this hmit is the same as the one for flat rectangular polyhedra. In particular this holds 

for the function S of flat wedge sets, what corresponds to n = 1 in the last equation. 

Now consider the construction of fig. (12) involving d — n + 1 dimensional polyhedra. 
The set Xi is formed by the union of two copies of a fiat rectangular polyhedron U = 
R X [a] having principal planes that form a small angle 9. Using the Lemma 2 we have 
that 7r^(^2) ^ 75"^(^i) + 2svol(i?)5j* in the limit of small a (see fig. (12)b). By 
the same reason 7^" ^(^2) ^ ^s~\^3) + 2svol(i?)(5^ in that limit ( see fig.(12)c). Then 
it follows that lima^o 7s'~^(-^i) = linia^o 75~^(-^3)- This means that we can move each 
of the components that are copies of U in Xi a direction perpendicular to their principal 



3. j{X) is proportional to the volume 




(67) 
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(a) 



(b) 



FIG. 11. Geometrical constructions used in the text to show that the hmit of the entropy for 
a flat wedge-like d — n + 1 dimensional set T coincides with jg{R), where R is the rectangular 
polyhedron forming the big phase of T. The picture shows a transversal view, where R is reduced to 
an interval, (a)- Two copies of T are placed such that their intersection and union are rectangular 
flat polyhedra with big phase R (neglecting sets small in more than one perpendicular direction), 
(b)- Two rectangular flat polyhedra Ri and R2 are placed such that both Ri — R2 and R2 — Ri 
are equal to T modulo sets that are small in more than one direction. 

plane without changing the limit of the entropy for a ^ 0. Thus, it also follows that 

lim,^o7r'(^4) = hm,^o7r'(^3). 

Now we show that these two last limits have different expressions in terms of the function 
ag when ^ — > 0. With the help of the auxihary sets drawn with dashed and dotted lines in 
the fig(12)c and (12)d, and applying to these constructions SSB and SSA respectively we 
obtain 

lim7r'(^3) < aU<^i,...,a,.n)+'2sSyol{R), (68) 
lim7r'(^4) > dg(2ai,a2,...,ad_„) +2s5«vol(i?), (69) 

a— >0 

where we have written explicitly the side lengths of the rectangular polyhedron that 
enter in the arguments of and (ai, arf_„) are the side lengths of R. Therefore 
Q;g(2ai,a2, ...,ad-n) = as{ai,a2, ...,ad-n), and from here 

Oisiai: aa-n) = a^{ai/2'', aa-n) , (70) 

for any integer r. In the limit of large r the right hand side of (70) is just Q;^^^(a2, dd-n), 
which is zero by the induction hypothesis. As the strongly subadditive an positive function 
a'g is zero for all rectangular polyhedra it must be zero for all polyhedra, since they can be 
written as union of simplexes, which in turn can be written as intersections of rectangular 
polyhedra. This completes the proof of the 

Proposition 2: Let S* be a li dimensional relativistic entropy function and a{X) = 
j{X) — 2s vol(X) the d — 1 dimensional Euclidean function defined in Section IV. If ao = 
then a = 0. 

4. Proof of theorem 3 

Finally, we prove the theorem 3. We take S with ag = and prove S{X) is proportional 
to the spatial border area of X. From proposition 2 it follows for the fiat polyhedra that 
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(a) 




X4 



(c) (d) 

FIG. 12. The entropies of Xi and are equal in the hmit of small transversal width of their 
connected components. This can be seen by constructing the auxiliary set X2, and applying Lemma 
2. More generally, we can translate any of the components of X3 in a direction perpendicular to 
its principal phase and get a set with the same limit entropy. This shows that X3 and X4 have 
the same limit entropy. The construction in (c) involves a set formed by two parallel rectangular 
components, the upper component shown with dashed lines and the lower one coinciding with the 
lower component of X3, and two wedge sets shown in dotted lines. The application of SSB leads to 
(68). The construction in (d) involves the set X4, two wedge sets and a flat rectangular set shown 
with dashed lines. Application of SSA leads to (69). 



S{X X [a]) — > 2s vol(X) for any d — 1 dimensional X in the limit a ^ 0. Thus, in this case 
the theorem 3 holds, since the boundary area of X x [a] is 2vol(X). 

A d dimensional simplex W can be written as W = {{xi, ...,Xd)/{xi, ...jXd) = riVi + 
... + ra+iVd+i, with > for i = 1, ...,d + 1 and ri + ... + r^+i = 1}, where Vi, ...,Vd+i 
are the vectors pointing to the simplex vertices. We are interested in a simplex in which 
Vd+i — u + w, where the vector u lies inside the simplex phase F with vertices vi, ...,Vd and 
w is a vector perpendicular to this phase. We arc going to find an upper bound for the 
entropy S{W) in the limit of small \w\. First decompose W as union of the d simplexes 
formed with vertices u and Vd+i, plus any subset of d — 1 vectors among the vectors Vi, ...,Vd. 
We decompose each of these simplexes as union of wedge sets with small width a, plus a 
number of small sets, as shown in the fig. (13). In the limit |w| — > the previous results 
show that the entropy of the wedge sets is proportional to the d — 1 dimensional volume of 
their big phase. The entropy of the smaller sets is bounded by some constant times the d — 1 
volume of their projection on the plane of the simplex F, since these sets can be formed by 
intersections of flat wedge and flat rectangular polyhedra. 

Taking the limit of small \w\ and a we then have 

\im S{W) < 2s vo\{F). (71) 

\w\-*0 

Consider now any simplex U with vertices ui, ...,Ud+i on the hyperplane of time t — 0, 
and let Fi,...,Fd+i be their faces. Call v the centre of the inscribed sphere, the unique 
sphere that is tangent to all the faces of U, and let r be its radius. The point x with spatial 
component v at time t = r is at null distance to all Fi, Fd+i- Call Gi, with i = 1, d + 1 
the null simplex determined by the phase Fi and x, and Hij, with i ^ j, to the d — 1 
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a 



FIG. 13. A three-dimensional flat simplex W with vertices vi,V2,v^ and V/i. The picture shows 
the projection on the plane of vi,V2,V2- The limit taken in the text is when approaches a point 
u on this plane. The decomposition shown for the lower simplex of vertices vi,V3,V4,u is in terms 
of wedge sets of width a, plus some small sets that can be formed by intersections of wedge and 
rectangular flat sets. Thus, in the limit a — and V4 ^ u we have the inequality (71) from 
subadditivity. 

dimensional simplex that is the common phase of Gj and Gj (see fig. (14) a). A point y with 
spatial component v at time t — r — e would be at small distance from the faces of U for 
e small. The point y and a phase Fj of U determine a flat simplex G'^ in the above sense, 
which converge to Gi in the limit e — >■ 0. The union of the simplexes G[ for i = 1, ...,d + 1 
form a Cauchy surface for U. We can slightly contract each of the simplexes G- in their 
own hyperplanes to form the simplexes G'^, which are are at positive distance a from each 
other (see fig.(14)b). We can further impose to the faces of G" and G" that are adjacent to 
each other to be parallel. Then, we add flat elements K^j formed by the Cartesian product 
of simplexes H^j, slightly smaller than the Hij, and a small interval a + e', in such a way 
to fill in the gaps between the different G'- and G". We can take all the Kij and the G'- 
to be commuting and the intersections of the K^j with G'l and Gj to be non empty sets of 
the form if^ times small intervals. The union of all the G", the K^j, and a number of sets 
that are small in more than one perpendicular direction simultaneously and have vanishing 
small entropy when taking limits, form a Cauchy surface for U. In the limits y ^ x and 
a — the entropy for the ffat simplexes G'- becomes 2s times their big phase volume, which 
coincides with vol(Fi). The entropy for Kij as well as the one for its intersection with G'- 
or Gj converges to 2svo\{Hij) in this hmit (taking also H'^j — > Hij and e' — > 0). Thus, the 
application of SSA to this geometrical construction leads to the inequality 

(d+l d+1 i-1 \ 

^vol(F,)-EEvol(//i,) . (72) 
i=l i=2 j=l J 

Now, as the hypersurfaces Gi are null and have zero expansion, we have that the volume of 
Fi is equal to the sum of the volumes of the Hij for all j different from i. Thus, 

d+l d+l i-1 

J2yol{Fi) = 2j2E^ol{H,,) . (73) 

i=l i=2 j=l 

From here and the fact that J2i=i vol(Fj) = area(f/) by definition, we have 

S{U) < s area([/) . (74) 
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The formula (74) allows us to obtain the corresponding inequality for any rclativistic 
polyhedron X. This can be decomposed into simplexes Ui having empty intersection or 
intersecting in adjacent faces Hij. We use the same procedure as above. Slightly contract 
the simplexes Ui and fill the gaps in the Cauchy surface for X with sets Kij formed by 
product of simplexes H^j (converging to Hij when taking the limit) and small intervals. 
The Kij are taken to intersect with the Ui and Uj in sets that arc the product of H-^ and 
small intervals. Then, in the application of SSA to the construction, all the internal faces 
determined by adjacent simplexes Ui and Uj do not contribute. This is because, according 
to (74), the contributions from the entropies of Ui and Uj to the term corresponding to 
the phase Hij add to 2svol(ify), the contribution of Kij converges to 2svol(ify), while on 
the other side of the inequality, the intersection of Kij with Ui and Uj gives 4svol(i/y). 
Therefore, only the external faces of the Ui (those that are also faces of X) contribute to 
the inequality, and we have 

S{X) < s area(X) (75) 

for all X. 

Now take a rectangular polyhedron R of sides (oi, a,^), and call its faces Pi with 
i = l,...,d {R has two faces equal to each Pi). Consider also the sets P/ = Pi x [3a] 
and P/' = Pi X [a, a, a], where [a, a, a] is a one dimensional set of two components of size a 
separated by a distance a. We have lima^o S{Pl) = 2s vol(Pj) and lima^o S{PI') = 4s vol(Pi) . 
We can form a big rectangular polyhedron U by union of translated and non intersecting 
copies of R, the adjacent ones separated by a distance a, and filling the gaps by copies of the 
P/ that intersect with the copies of R in sets that are translated copies of P/'. To complete 
U we have to add sets that are small in more than one direction. Taking into account (75) 
the limit of S{U)/yol{U) for U going to infinity is zero. Thus, the application of SSA to the 
construction leads in the limit of small a and infinite U to 

S{R) > s area(P) . (76) 

Thus, the equality S{R) = s area(P) holds for any rectangular polyhedron. 

Any rclativistic polyhedron X is included into some rectangular polyhedron R. Complete 
a Cauchy surface for X with simplexes to form a Cauchy surface for R. Then, repeating the 
above procedure for dealing with the internal adjacent boundaries we obtain from SSA 

S{X) > s area(X) . (77) 

This finishes the proof of theorem 3. 

VI. DISCUSSION AND CONCLUSIONS 

The properties of the quantum entropy make possible to deal with the geometric entropy 
in a geometrical way. However, even if it is believed that this quantity has a deep physi- 
cal significance, it remains an abstract concept, wanting for a solid mathematical basis to 
compute it. In this sense, the principles employed in the paper to obtain the area law are 
very basic and general, and one can expect them to be still available once a finite geometric 
entropy is defined. 
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Hl2 

(a) (b) 

FIG. 14. (a)- The three dimensional set formed by the future domain of dependence of a two 
dimensional simplex U. The sets Fi, F2 and Fz are the faces of and the point x is at null 
distance from these faces. Gi, G2, and G3 are the faces of the null boundary of the future domain 

of dependence, and Hij is the common phase between d and Gj. (b)- A view of the set in (a) 
from the top. The simplexes G'l approach the null surfaces Gi and thus they are flat simplexes. 
With the simplexes G'l and the commuting rectangular sets Kij plus some small sets of negligible 
entropy limit, we form a Cauchy surface for U. 



A more immediate meaning can also be given to our results. They indicate that the 
entropy must be proportional to the area in the limit of infinite cutoff for any regularization 
scheme in QFT where covariant results are obtained asymptotically. In other words, the 
leading divergence in the entropy must be proportional to the area for any set. We remark 
that the theorem shown in this paper tells that no subleading term exists in the continuum, 
excepting the constant. 

This is indeed what all previous QFT calculations have shown for particular cases. Be- 
sides, we have also checked the area law for an spherical corona using the numerical method 
of ref. [5]. This case has also been studied in [19] in a related context. This is a very inter- 
esting example since it allows to test a very important and counterintuitive characteristic of 
the area solution, namely, that it is not ordered by inclusion outside the class of convex sets. 
As expected, the leading term in the entropy grows with the sum (with equal coefficients) 
of the areas of the internal and external surfaces of the corona. 

There is a curiosity about the Euclidean solutions which could be related both with the 
uniqueness of the relativistic entropy function and to the divergences. In the Euclidean case 
it is possible to construct many entropy functions averaging over translations and rotations 
positive and strongly subadditive functions which are not Euclidean symmetric. Among all 
these, we can find in particular some very interesting monotonic solutions that arc exactly 
proportional to the area for the convex sets and less than the area for the non convex ones, 
and where the entropy of two distant objects 'interact' by their mutual shadow. We suspect 
that a classification of Euclidean solutions could be possible using this averaging method. 
Coming back to the relativistic case we see that a similar construction does not work since 
the Lorentz group is non compact, leading to divergent integrals. This is so even taking 
compact sets to make finite the integration over translations. 

An interesting question related to the QFT calculations is the physical meaning of the 
coefficient relating area and entropy and the constant term. The coefficient of the area term 
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was related to the Newton constant in [6,23] (see however [8,10,24]), making contact with 
the black hole entropy. For 1 + 1 dimensional conformal field theories it was shown in [9] 
that the most divergent term for the single interval entropy is a constant (independent of 
size) proportional to the central charge. This coefficient corresponds to the combination 
7 = a + /3 according to (54). A natural question here is if a new parameter shows up 
in the expression for the multicomponent set entropy [25]. In other words, the entropy 
can be simply proportional to the number of components {a = 0), or it can contain an 
additional constant term a, with the same degree of divergence as /3. It is possible that this 
last term is always subleading and negligible in the continuum. This is specially expected 
in more dimensions since the constant term requires a dimensionless parameter, while the 
dimensionfuU one in front of the area term diverges with the cutoff. Thus, the constant may 
represent a spurious solution for c? > 1. 

The entropy per unit area s is a function Sp of the Poincare symmetric state p. The 
combination of the inequalities (2) and Theorem 3 imply that given any two such states we 
have 

= AiSpi + A2Sp2 , (78) 

for positive Ai and A2 with Ai + A2 = 1. Thus the entropy per unit area is afin under mixing 
of states, a property known for the mean entropy per unit volume of traslational invariant 
states. The physical meaning of eq. (78) is however unclear to the author. 

In a series of recent papers a very interesting phenomena has been pointed out [11,26]. 
The authors discovered that the entropy in a volume is not the only quantity that scales 
with the area of the given region but the same behavior is expected for a whole class of 
operators averaged using the local density matrices. In particular they showed the area 
law for the energy fiuctuations. The results of this paper rely heavily on the use of the 
strong subadditive property of the entropy, and thus they can not be extended directly to 
operator correlations. However, it is possible that the method here could be generahzed by 
introducing chemical potentials for suitable operators attached to the local regions. 

The proof in this paper uses explicitly arbitrarily small and large distances. The infinite- 
ness of the space is specially relevant in the argument around the constructions of fig. (12), 
which involve large translations, and lead to the equation ^{X) = (3{X). This is what ex- 
cludes the possibility of finding a solution ordered by inclusion. However, it is not excluded 
that different type of solutions exist in de Sitter space. This has compact Cauchy surfaces, 
and it is also maximally symmetric, what most probably can allow us to make a complete 
classification of solutions as the one given here for the Minkowski space [25]. 

The entropy of the black hole and the geometric entropy are 'empty' space quantities. 
There are a number of proposals for general bounds to the entropy of bounded systems 
which are suggested by the physics of black holes. They have been successfully checked 
mainly in the classical regime. The definition of a bounded system requires a partition of 
the Hilbert space into a tensor product, and it is necessary to add boundary conditions. 
Thus, at the semiclassical level, the entropy bounds should include a term coming from the 
vacuum correlations, even in empty space (see similar ideas in [20]). In principle, this could 
greatly tighten the inequalities. How do these bounds compare with our result in flat space?. 
Consider the covariant bound proposed in [27], which is related both with the Bekenstein 
bound [28] and the generalized second law. Take a d — 1 dimensional spatial surface fli 
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and a mill congruence W orthogonal to it, having non positive expansion, and ending in 
another ri — 1 dimensional spatial surface The bound states that the entropy crossing W 
is bounded above by a constant times vol(r2i) — vol(r22). Thus, this gives zero entropy for 
our polyhedra, since they have Cauchy surfaces formed by union of null fiat surfaces having 
zero expansion. 

This touches upon the question of the geometric entropy for a more general type of 
subsets than the one considered here. The area and the constant term are not the most 
general solutions of the SSA inequalities for larger classes of subsets than the polyhedra. 
Take for example the class of subsets with piecewise differentiable spatial corner. It is also 
a closed class under intersection and union of commuting sets. As can be easily seen from 
fig(l)b, the integral over the spatial corner of any geometrical quantity that is positive and 
local gives place to a relativistic entropy function. Of course, the theorem 3 implies that 
these solutions are proportional to the area on the polyhedra. As an example of such a local 
quantity we can take the square of the expansion of the future null boundary of the given 
causally closed set. In this case the entropy is non zero for the sphere and zero for every 
polyhedron approaching it. This shows that this type of solutions has a very pathological 
discontinuous behavior. 

Excluding these solutions, probably unphysical, it seems we are left with two possibilities, 
either the covariant bound is wrong outside the classical domain, or the assumptions of this 
paper fail somewhere. However, it is also possible and likely that the covariant bound has 
to be applied to a regularized notion of entropy which is zero on Minkowski space. This is 
appealing, since the calculation of the black hole entropy using the Euclidean path integral 
requires the subtraction of the corresponding flat space quantity [29] . 
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